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The connection between quantum optical nonclassicality and the violation of Bell's inequalities 
is explored. Bell type inequalities for the electromagnetic field are formulated for general states of 
quantised radiation and their violation is connected to other nonclassical properties of the field. This 
qI^ , is achieved by considering states with an arbitrary number of photons and carefully identifying the 

■ hermitian operators whose expectation values do not admit any local hidden variable description. 

We relate the violation of these multi-photon inequalities to properties of phase space distribution 
functions such as the diagonal coherent state distribution function and the Wigner function. Finally, 
the family of 4-mode states with Gaussian Wigner distributions is analysed, bringing out in this 
case the connection of violation of Bell type inequalities with the nonclassical property of squeezing. 
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I. INTRODUCTION 



> 

\ The superposition principle and the reduction of the wave function upon measurement are the two essential features 
of quantum mechanics, which are responsible for most counterintuitive paradoxical situations that arise when one starts 
to interpret the predictions of quantum mechanics. The violation of Bell's inequahties is one of quantum theory's most 
qI^ : striking consequences Q. When outcomes of individual measurements are given an objective meaning, the locality 
condition of "no action at a distance" imposes constraints on the possible correlations, which are expressed through 
Bell's inequalities. There exist quantum mechanical states for which these inequalities are violated bringing out the fact 
that in quantum mechanics we either have to give up hope of an objective interpretation of individual measurements 
or accept "action at a distance" . The possibility of a local "classical" theory of individual measurements is thus ruled 
^ out. Though initial work, starting with Bohm concerned itself with the states of two spin-i particles, it was for 
^ ' particular states of the quantised electromagnetic field that experiments were first performed in this direction Q [|j . 

The formulation of Bell-type inequalities is possible only for a system which has two or more kinematically indepen- 
^ ' dent subsystems. These subsystems could each, for example, be a spin-half system with states in a two dimensional 
• i-H ■ Hilbert space, a photon with fixed energy-momentum but variable polarisation state, or even a quantum mechanical 
' system with one canonically conjugate pair of operators q and p and states in an infinite dimensional Hilbert space |^ 
@ 0- Many discussions focus on states of two photons, for two fixed propagation vectors and variable polarization 
- - ■ states. To go beyond a single photon in each mode, in our analysis we will deal fully with all states of 4-mode fields, 
which include the two-photon states as a simple special case. On the other hand, the two-mode electromagnetic fields, 
with different directions of propagation and fixed polarisations, provide us with a situation kinematically equivalent 
to the one considered in the EPR paper as the subsystems are of the q,p form. 

The quantum mechanical states for which Bell's inequalities are violated thus have essential "quantum" features and 
defy a deeper interpretation based on realism (the objective existence of attributes independent of their measurement) 
and locality (no action at a distance). On the other hand, for a state which obeys a complete set of Bell's inequalities 
it is in principle possible to give a local "classical" interpretation to individual measurements. An independent notion 
of classicality, used in quantum optics ||^, for the states of the quantised electromagnetic field is based on the 
diagonal coherent state distribution function, and here, the specific classical theory one has in mind is Maxwell's 
theory. This paper explores the connections between these two different ways of classifying the quantum mechanical 
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states of the electromagnetic field as "nonclassical" and "classical". To carry out such an analysis, we obviously 
have to go beyond two-photon states and formulate Bell type inequalities for general four mode states with arbitrary 
photon number distribution. 

A typical setup used to study violation of Bell's inequalities involves four modes of the field with propagation 
in two different directions, and arbitrary polarisations being allowed transverse to each direction. For photons in 
each propagation direction a particular polarisation is selected by a variable polariser, and finally coincidence counts 
are recorded using photo detectors. In order to analyse an arbitrary state of the 4- mode field and derive a Bell- 
type inequality for it, two crucial inputs are needed: the first is the identification of hermitian operators whose 
measurements and correlations are directly related to the coincidence count rates; and the second is an appropriate 
quantum mechanical analysis of the polariser, namely its action on a general input multiphoton state. We feel that 
previous analyses in this direction lack proper treatment of one or the other of these aspects Q| 

A coincidence count may be said to be registered when one or more photons are detected at each of the two 
detectors simultaneously or within a pre-assigned time interval, disregarding the exact numbers of photons detected. 
(So this is a kind of coarse grained coincidence rate). Hence such coincidence count rates differ from intensity- intensity 
correlations. We introduce an appropriate set of hermitian operators in the space of states of the 4-mode field so that 
all such coincidence counts can be calculated from the expectation values of products of these operators. By assuming 
a local realist hidden variable description for the outcome of various eigen values of the above mentioned operators 
we derive Bell type inequalities constraining the measurable coincidence count rates. 

Classical states in the quantum optical sense have an underlying classical distribution function; they have been 
shown not to violate Bell's inequalities, while nonclassical states are potential candidates for such violations. The 
quantum optical nonclassicality is invariant under passive canonical transformations which form the group U (n) for 
a general n mode system p3| ; such transformations can be experimentally implemented in a straightforward way 
using optically passive elements like beam splitters and mirrors. On the other hand, such transformations can alter 
the potential of a state to violate Bell type inequalities. Therefore, starting with a nonclassical state we should allow 
it to undergo arbitrary passive canonical transformations before looking for violations of Bell's inequalities. It may 
turn out that a state which obeys Bell's inequalities is related by a passive transformation to one which violates them; 
such a violation obviously is a consequence of the nonclassicality of the original state. This is closely connected to the 
fact that passive transformations can take nonclassical nonentangled states to nonclassical entangled ones. Generally, 
therefore, the relationship between total photon number conserving passive canonical transformations on the one 
hand, and entanglement or violation of Bell-type inequalities on the other, has to be carefully analysed. We show that 
this capacity of U (n) to alter the entanglement properties of states underlies the violation of Bell type inequalities by 
beams originating from independent sources discussed by Yurke and Stoler Further we show that for their 

scheme to work at least one of the beams has to be in a quantum optically nonclassical state. A much larger set 
of quantum mechanical states of the 4-mode electromagnetic field can be analysed in our formalism opening up the 
possibility of more easily experimentally observing the predicted violations of Bell's inequalities. 

The material in this paper is arranged as follows: In Section II we describe the basic 4-mode setup to be used for 
the study of Bell's inequalities. The coincidence count rates for a general 4-mode state are then defined and used 
to derive Bell type inequalities with an appropriate quantum mechanical description for the polariser. Section III 
explores applications to various 4-niode states. Coherent states which are "classical" in more than one way {i.e., they 
have classical diagonal coherent state distribution functions, are not entangled in any basis and lead to classically 
expected "photon counting"), are shown to obey Bell type inequalities. This result is then used to prove that all 
quantum optically classical states also obey our Bell-type inequalities. Finally we analyse in detail the family of states 
of the 4-mode field with a Gaussian- Wigner distribution, bringing out their violation of the Bell-type inequalities and 
the connection of such violation with the nonclassical property of squeezing. Section IV contains some concluding 
remarks. 



II. BELL'S INEQUALITIES FOR MULTIPHOTON FIELDS: CHOICE AND PREPARATION OF STATES 

Consider four modes of the electromagnetic field chosen so that we have two different propagation directions labeled 
by wave vectors k and k' and along each direction we allow arbitrary polarisations. We can choose for convenience 
a basis for each polarisation space and can then label these modes with annihilation operators 01,02,03 and 04; the 
modes ai and 02 refer to the linear polarisations along x and y for the beam represented by wave vector k and the 
modes 03 and 04 to the linear polarisation modes along x' and y' for the wave vector k'. Without any loss of generality, 
we can assume that both the directions of propagation are in the plane of the paper; then x, x' are also chosen to be 
in the same plane while y, y' are in the common direction pointing out of the plane of the paper. These are indicated 
in the left end of Figure 1. The t/(4) block symbolically represents the replacement of oi, 02, 03, 04 by complex linear 
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combinations of themselves determined by a matrix of J7(4): 

a'j = Y.Ui''''k, UgU{A) (2.1) 

k 

which is a passive total photon number conserving canonical transformation. Polarisers Pi and P2 can be set at any 
angles 0i and 02 with respect to x and x' respectively. After passing through these polarisers the beams encounter 
detectors Di and D2 connected to each other by a coincidence counter. This scheme is very similar in spirit to most 
experimental and theoretical situations which have been extensively investigated. 




C.C. 



FIG. 1. Setup to study the violation of Bell type inequalities for arbitrary states of the four mode radiation field. The U(4) 
block represents a variable passive U (4) canonical transformation mixing the four modes, Pi , P2 are the polarisers oriented at 
angles 61,62 with respect to the axes x and x' respectively. Di and D2 represent the detectors and the block c.c. indicates the 
coincidence counter. The part of the diagram to the left of the vertical line indicates symbolically the state preparation stage, 
that to the right the measurements done on it. 

However, the difference here is that we do not wish to assume anything about the incident 4-mode state; it could 
have arbitrary number of photons, and could even be a mixed state. The measurements of interest here are the 
coincidence count rates, either for given settings of the polarisers {i.e.,9i and 62) or when one or both polarisers 
are removed. For a given 0i{02), Di{D2) receives and detects photons in the corresponding linear polarisation state 
alone. (For simplicity alone we restrict the analysis here to linear polarisation states) 

A. Interpretation of the coincidence count rates 

A coincidence is counted wlien both the detectors Di and D2 click simultaneously('Simultaneous' here has the 
meaning as stated in the introduction that within a preassigned time interval At each of the two detectors registers 
at least one photon). We will not distinguish in this analysis between coincidences of different strengths i.e., different 
numbers of photons (greater than one) being received by each detector. Coincidence rates defined in this specific 
physical manner must be represented by suitable choices of hcrmitian operators and their expectation values; this will 
be done in detail below. In our analysis, we will need the following four types of coincidence count rates; 

(a) P{0i,02)- The first polariser at 9i and the second one at 02 with respect to their respective x axes. 

(b) P{0i, ): The first polariser at 0i and the second one removed. 

(c) P{ , 62): The first polariser removed and the second one at 02. 

(d) P( , ): Both the polarisers removed from the setup. 
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These coincidence counts are the measurable quantities and can be calculated for any quantum mechanical state of 
the 4-mode field. On the other hand, we will see that these count rates have to obey certain inequalities if we demand 
a realist description and locality as well, so that measurements at detector Di are independent of the setting of the 
polariser P2 and vice versa. 

B. Quantum mechanical description of the polariser 

Classically, the action of a polariser is straightforward. It allows a particular polarisation state (which defines 
the polariser axis) to pass through without any hindrance and blocks the orthogonal one completely. For a general 
polarisation state, the component of the electric field along the axis passes through unaffected while the orthogonal 
component is absorbed. We are interested in the generalisation of this action to quantum mechanical situations. For 
a comparable discussion for the beam splitter action see [^. Given the classical action, we know how the polariser 
should act on single photons and on coherent states. For a single photon traveling in the z direction, in an arbitrary 
state of polarisation ci|l)a;|0)j, + C2|0)a:|l)y (ci, C2 being arbitrary complex constants with |cip + |c2p = 1)0 in the x — y 
plane and with a polariser placed in the x direction, the probability of transmission is |cip. From this we can build 
up the probabilities of transmission for arbitrary states as they can be expanded in terms of number states. Consider 
a general pure input state n2 ^nin2\'ni)x\n2)y with the polariser placed along the x direction. The probability 
of finding ni photons after passage through the polariser should be given by |C'nin2p- Nevertheless, this still 
leaves the question about the output state open-ended; what is the state of the single-mode field which emerges 
out of the polariser? Exactly identical probabilities can be obtained from pure as well as mixed states emerging 
from the polariser. Compared to the classical case, we need to give more detailed consideration to the physical 
process underlying the removal of photons of particular polarisation from the beam; it could either be absorption by 
certain other degrees of freedom(lattice) or a change of direction caused by different refractive indices for different 
polarisations. In either case, after passage through the polariser the information contained in that mode, though 
existing, is inaccessible and therefore we should trace over that mode to obtain the outgoing state of the field. 

We thus arrive at the following general prescription for the action of the polariser: for a given input two-mode state 
with density matrix p, the two modes being two orthogonal polarisations along the same direction of propagation, 
and a polariser placed at an angle 6 with respect to the x axis, the single-mode state p(9) after passage through the 
polariser is obtained by taking the trace over the mode orthogonal to the linear polarisation defined by 9: 

00 

Here we have chosen the number state basis for the mode orthogonal to 6; we could as well have chosen any other 
complete set of states. More explicitly, consider the most general density matrix for the two polarisation modes along 
the same direction of propagation. 

P= ^ Cn^n2n[n'Jni)x\n2)y x{n'i\y{n'2\ (2.3) 

If a polariser is placed at an angle 9 with respect to the x axis we first rotate the basis in the x ~ y plane to 9 and 
9 + ^ giving us 

P= X! C'«in2n;^(6')l»^i)el"2)e+f ei^^ile+f ("2! (2.4) 

Each C„j„2„'^„^ (0) is a dependent linear expression in the C....'s, of course conserving ni + n2 and n[ + n'2. Now 
tracing over the mode orthogonal to 9 yields the final state after the passage through the polariser: 

P(^^^ E (E^"i"^S«2(^)) Moo{n[\ (2.5) 

ni ,n'^ \ '^2 / 



^Here we have used the language appropriate to two modes and allowed the single photon to be in either of the polarisation 
modes; equivalently we could have talked about a single photon in a superposition of two orthogonal polarisation states 
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When p is a pure state «.e., C„^„2„'^„^ = Cnmj^n'n'i the rotated coefBcients C„j„2„'^„^(^) also factorise, and the 
output state 

/'W = E (E^"i"=W^nin.w) \ni)ee{n',\ (2.6) 
could in general be mixed. For the special case of coherent states |2;i)x|-22)j/ 

^ _ -^(i^ii' + i^2r) 4^ 



n -^(kil' + l^2n (2icos0-z2sin0)"i (2isin0 + Z2COS0)"^ ,„ 

^nin2\y)=^ 1=. 1=^ V^-'J 



ni! \/n2! 



Therefore the final density matrix after the beam emerges from the polariser is given by 

Zl\^ + \Z2\'') V- (£lC0S^^£2Sin^ 



(4cos6'- ^sin^)"! v-^ (|0isin6i + 02cos6'|)^"= , 



_,-|zicos^-Z2sin^|2 ^ (zicosg-02sin^)"^ {z^cosO - z^sm9)< , 

ni,n5^ V i VI 



(2.8) 



which is a density matrix corresponding to the pure coherent state \zi cos 9 — Z2 sin9)g, consisting of photons all of 
which are polarized in the direction 9. The reason for this is the subtle fact that coherent states are not entangled 
in any basis we may choose for the two-mode incident state, and thus tracing over one of the modes just amounts to 
neglecting the corresponding factor in the product state after an appropriate rotation of the basis 

|-2i)x|-22)y = l^^i cos^ — ^2 sin^)^!^! sin^ + ^2 cos^)e+:| — > cos^ — cos^)^ (2.9) 

On the other hand, single photon states can in general be entangled states of the two-mode field and thus would 
lead to mixed one-mode states after passing through the polariser. For example a pure two-mode single photon state 
-^(|l)j;|0)j, + |0)j;|l)j,), after passage through a polariser placed in the x direction reduces to a mixed state with 

density matrix ^{\0):ox{0\ + |l)xx(l|)- 



C. Derivation of the inequalities 

In order to define and calculate the coincidence count rates, consider the following four hermitian operators, all 
having eigen values and 1 

^l=(/2x2-|00)(00|)k 
l2 = (/2x2-|00)(00|)k, 

A,{9,)^{Ie,-\0)e, e,(0|)/e,+f 

M02) = {le, - \0)e, eM) Ib,+^ (2-10) 

The subscripts 9i and ^2 in the last two equations refer to the directions of the polarisers. Thus Ai and Ai{9i) are 
operators belonging to the first two modes of our 4-mode system, namely propagation direction k and polarisation 
along X or y; 72x2 is the unit operator for these two modes while le-^ is the unit operator for the single mode propagating 
in direction k and with polarisation 0i, and /e2+f that for the orthogonal polarisation. A2 and ^2(^2) are defined in 
a similar way for propagation direction k'. Expectation values of the above operators are probabilities of finding at 
least one photon of the appropriate kind: 
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(^l(^l)) 
{MO2)) 



probability of detecting at least one photon at Di with Pi removed, 
probability of detecting at least one photon at D2 with P2 removed, 
probability of detecting at least one photon at Di with Pi set at Oi, 
probability of detecting at least one photon at D2 with P2 set at ^2- 



(2.11) 



These operators are defined so that when we calculate their expectation values in any state, the actions of the 
polarisers as described in the previous subsection are automatically implemented. In an individual measurement, 
one of the eigen values of the hermitian operator is recorded, with the probabilities for the two possible eigen values 
being calculable from the wavefunction or the density matrix. On the other hand if a hidden variable description is 
available, then given the value of the hidden variable A as well as the quantum mechanical state vector, we can in 
principle predict the outcomes of individual measurements. It is well known that while keeping the hidden variable 
theory completely general and by imposing the locality condition (no action at a distance), general constraints on 
the correlations between observed quantities can be obtained; we will perform a similar analysis for our situation and 
derive Bell type inequalities for the correlations among the A's. 

Assuming the probability distribution for A to be P(A)(for simplicity the dependence of V{X) on the specific 4-mode 
quantum state is omitted) so that 



we now proceed to develop expressions for the average values of various observables. Let ai(A),a2(A),ai(^i, A), and 

02(^2, A) be the actual values of the dynamical variables Ai, ^2, v4i(0i) and A2{62) for a given value of A respectively: 
they take values or 1. The averages computed from hidden variable theory are then given by 



The average values of the products of these operators, corresponding to the various coincidence count rates, are then 
given in local hidden variable theory by 



The subscript "hv" is for 'hidden variable' and distinguishes these averages from the quantum mechanical expectation 
values. Here in the first step, ai2's are the hidden variable theory values of corresponding dynamical variables A1A2 
etc. The last step in each of the above equations is very crucial and is basically the expression of the locality assumption 
on the hidden variable theory. We have assumed for instance that the quantity 012(^^1, f?2, A) is a product of the two 
factors ai{di, A) and 02(^2, A), each depending only on one angle and not the other; more generally what we measure 
for the propagation direction k does not depend on what we choose to measure (or not measure) for direction k'. This 
assumption is "reasonable" because we can arrange the situation such that the two measurement events are space like 
separated. 




(2.12) 




(2.13) 




(2.14) 
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One more assumption we make is the "no enhancement assumption"^ i.e., the presence of the polariser can only 
remove photons from the beam and is incapable of adding photons or increasing the coincidence count rate: 



ai(0i,A) < ai(A) 
02(^2, A) < a2(A) 



(2.15) 



We now recall the lemma due to Clauser and Horne which we will use to derive inequalities for the above calculated 

correlation functions. 

lemma: 

if < X, a;' < X and < y, ?/' < r then 

- XY <xy- xyf + x'y + x'y' -Yx' - Xy <0 (2.16) 



With the "no enhancement assumption" (2.15) and the product forms (2.14), after identifying x = ai(0i,A), x' — 
ai{9[, X),y — 02(^2, A), y' = 02(^2; A), X = ai(A), Y = 02 (A) we get the following inequality 

- ai (A)a2 (A) < ai (O^ , A)a2 (^2 , A) - ai (^i , A)a2 (^?2 , A) + ai ie[ , A)a2 (^2 , A) + ai ie[ , A)a2 {02 , A) 

-a2{X)ai{9[,X) - ai{X)a2{02, A) < (2.17) 

Integrating the above inequality over A, with weight function V{X) we get the inequality obeyed by the coincidence 
count rates for any choices of angles 6*1 , 6*2 , O'l , 6*2 

, )<Piei,e2)-P{ei,e'^) + p{e[,e2) + P{e[,e'^)-p{0[, )-p( ,62) <o (2.18) 

This is the inequality we will use in our future analysis. The left hand side inequality depends upon the total 
coincidence count rate without the polarisers whereas the one on the right hand side does not. This result is of the 
same form as that derived and used by all workers on the subject; what we emphasize here is that we have clearly 
specified the hermitian operator observables at the quantum mechanical level, for which correlations are considered, 
so the actual expressions for various P's are specific to our treatment. 

For such a theory to be consistent with quantum mechanics, the coincidence rates computed from a quantum 
mechanical calculation should also obey the above inequalities. In the above derivation we have not assumed anything 
about the 4-mode state and therefore we can consider arbitrary states to check for any violation. 

If a given quantum mechanical state does not obey the above inequalities then it definitely has some nontrivial quan- 
tum features which cannot be accommodated within realist hidden variable theories based on locality. We now discuss 
various examples of 4-mode quantum states in sequence, and examine the validity or otherwise of inequalities (2. IS) 
in each case. 



III. MULTIPHOTON STATES, VIOLATION OF BELL'S INEQUALITIES AND PHASE SPACE 

DISTRIBUTIONS 

A. Two-photon states 

Before undertaking the analysis of multiphoton states, we consider here the extensively studied two photon state 
as a warm up exercise. It so happens that for these states our formalism reduces to the usual one and we get results 
identical to those already available in the literature. Consider the following pure two-photon state of the 4-mode field 

W = ^(|l)l|0)2|0)3|l)4 - |l)l|l)2|0)3|0)4 - |0)i|0)2|l)3|l)4 + |0) 1 1 1) 2 1 1) 3 1 0) 4) 

= i(al - 4)(4 - 4)|0)i|0)2|0)3|0)4 (3.1) 

The calculations of the quantum mechanical coincidence count rates become easy if we make the following obser- 
vation about the operators Ai{9i) and ^2(^*2): 



^This is at the level of hidden variable description. We have already seen that the possible values of ai{di, A) and ai are or 
1. So if ai(A) = for some A then ai{9i, A) = as well for any 9i. 
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A2{e2)=u{R2{e2)) ih - \o)3 3{o\)hu-\R2{e2)) (3.2) 

Here U{Ri{9i)) and W(i?2(02)) are the unitary operators corresponding to the transformations Ri{9i) and i?2(^2) 
which rotate the basis in the polarisation space by angles 9i and 82 for directions k and k' respectively. 

R^{0,)=( '""'^ RM = ( '""'^ ''''I') (3.3) 

^ ' smt'i coswi J ^ ' sin^2 cost^2 / 

Using this form of the operators the calculation of the expectation values of various combinations of operators is 
straight forward and yields the quantum mechanical coincidence count rates: 

i sin^ {B^ + 02) 





^'2)q?Ti — 




P{di, 


stp 


{m,{0,)A2m = \ 


Pi , 




{i,\MA2{02m = \ 


P{ , 


stp 

/qm ^ 





(3.4) 



(The superscript 'tp' here refers to two-photon states). Substituting these rates in the inequality ( 2.18 ) we get the 
following condition, if this state is capable of a local realist hidden variable description: 

- 1 < i(sin2(6'i + 62) ~ sin2(e'i + e'2) + s\n^{9[ + 62) + s\n^{e[ + S^) - 2) < (3.5) 

We find that the above condition can be violated on either side for some values of 9i,9[,92 and 92- As an example 
let us choose 9i = ^,92 — ^,9[ = ^ and 9'2 = 0; then the right hand side of the above inequality becomes 

^{V2-l)<0 (3.6) 

which is clearly violated. 

Since two-photon states have been extensively studied in the literature we connect our result with existing analyses 
and make some useful observations. The single hermitian operator A = / — |0)(0| whose expectation value gives the 
probability for finding one or more photons reduces effectively to a' a for the appropriate mode when at most one 
photon is present in the beam; therefore our results agree with existing ones For two-photon states we have the 
relations P{9i, ) = P{9i, 6*2) -I- P{9i, 6*2 + f) etc. They too are a consequence of the reduction of the A operators to 
the above described form and do not remain valid for the analysis of general states. 

B. The coherent states 

As a first nontrivial example we consider 4-mode coherent states defined as (omitting further subscripts outside 
kets) 

|zi)|z2)|z3)|z4)= e ^ V y X 

eV ' 1 '2 33 V 10000) (3.7) 

where zi, Z2, 23, Z4 are complex numbers. The calculation of quantum mechanical coincidence count rates for this case 
is rather straightforward and we get the following results (the superscript 'cs' means coherent states): 

Pi9,,92)^,n = {l-e-\'i\"){l-e-\4\") 

Pi&i, )qm = (1 - e-|^il')(l - e'^^'^^l^ + I^^H) 
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P{ 


> ^2)qm 


= (l-e 








\CS 


= (1-e 








^1 


= Zi cos 


7i — Z2 sin c/]^ 






4 


= zi sin ^ 


'1 + Z2 COS 6'i 








= Z3 cos 


92 — Z4 sin 6*2 






4 


= Z3 sin ^ 


h + Z4 cos 6*2 





(3.8) 

We see that P^^-^{di, 62) is a product of two factors, the first one depending on 61 alone and the second on 62 alone. This 
fact originates from the unent angl ed nature of the coherent states and will be used now to show that coherent states 



always obey the inequalities ( 2.18 ). Using the facts that < P(9i 
P{ 1 )qfn and the lemma ( 2.16 ), with the identification x — P{9i 
P{ , 6'2)q|n' X = Y = P{ , '^^ arrive at the following preliminary inequality 



\cs 
^qm 

NCS 



cs 



x'^P{e[, )cs 



qmi Pi 
qmi 



2Jam) Pi 



qm 

y = P{ 



Ql \CS < 

-'2^qm ^ 



92)qmi 



'P{ 



\2 CS ^ 

jqm — 

|CS p/ 

qm-* V 



P{(^i: Jqm 

+P{0'i. 
< 



lqm^( 



i^2)qm -P(^i: )qrn^( i^2)qrn 
' ^2)qrn ^ )qrn-P( 



';)qrn + ^(^ii )qrn-P( : 
)qrn ^ P{ i ^2)qrn^'( 



^2)qm 

\CS 



Dividing throughout by the positive number P( 
initial inequality reduces to the Bell form: 



)qm ^^"i using the expressions of eq. (|3 



(3.9) 

I) we see that the above 



Pi I )qm < -P(^'i)^2)qm --P(^i)^2)qni + -P(^'ii^2)qm 



^cs 



qm 



Pi ,^2)^?n<0 



(3.10) 



Thus we s ee ex plicitly that coherent states do not violate the inequalities (2.18). The reduction of equation (3.9) to 
equation ( ^3.10 ) is possible because coincidence count rat es for coherent states factorize. Further, we emphasize that 
it is p ossib le to write an inequality corresponding to ( |3.9| ) for any quantum state but it can be reduced to a Bell type 
form ( 3.10| ) only for those special states for which coincidence count rates factorize. 



C. The classical states 

We can express any arbitrary state of the 4-mode radiation field in terms of projections onto coherent states in 
particular, for a state with density matrix p, we have 

P = ^ I Vizi, Z2, Z3, Zi)\zi, Z2, Z3, Z4){zi, Z2, Z3, Z4\(f Zl(P Z2(fz3tPzi, (3.11a) 

■\ J ip{zi, Z2., Zz, Z4)<P ZlCp Z2(P Z'id? Zi — 1 (3.11b) 
TT 

The diagonal coherent state distribution function (^(zi, Z2, 23, 24) describes the state p. In quantum optics the states 
of the 4-mode field are classified into classical and nonclassical types as follows: A given state is classical if the 
diagonal coherent state distribution function for it is nonnegative and nowhere more singular than a delta function. 
Otherwise the state is nonclassical. The function undergoes a point transformation when the state undergoes a 
unitary evolution corresponding to a passive canonical tranformation given by an element of U{A): 

V3(Z1, 22,23,24) p izi,Z2,Zj„Z4) = <y9(zi,Z2,4'4) 

{z[ z'2 z'^ z'^) = (zi Z2 Z3 Zi)U^ , U(^U{A). (3.12) 

Therefore, the classical or nonclassical nature of a state is preserved under such transformations. 

Coincidence count rates, and for that matter expectation values of any observable can be calculated from the 



diagonal coherent state distribution function. For a classical state, the inequality (2.18) becomes, by multiplication 



of all terms in (3.1C) by p followed by integration 
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)qin'i5(^i: 2^2, 23, Z4^)(f zi(f Z2(f Z3(f Z4 < 



7? y"'-^*^^i'^2''qm--P(^i7^2)qni + -P(^i!^'2)qm+ )qm 



-P( , 02)'i(^in)'Pi^li ^2, Z3, Zi)(f Zitf Z2(f zscf Z4 < 



(3.13) 



The above result makes use of the nonnegative nature of ip and will not be true for a nonclassical state. Using the 
normalisation (3.11b), and the nonnegative nonsingul ar na ture of 93(^1, Z2, Z3, 2:4), we conclude that a "classical state" 
will not violate Bell type inequalities defined in eq. ( 2.18 ). Further, since the classical or nonclassical status of the 
4- mode state is invariant under passive canonical transformations which form the group C/(4), a classical state after 
undergoing such transformations will still not violate Bell type inequalities. On the other hand, the nonclassical states 
i.e., the states with negative or singular diagonal coherent state distribution functions can violate these inequalities; 
in fact, the violation of such an inequality implies that the underlying diagonal coherent state distribution function 
for the state is negative or singular i.e., the state is nonclassical in the quantum optical sense. 



D. Application to two-mode Gaussian states 



We have seen that typically a minimum of four modes are required for the analysis of Bell's inequalities; but 
if we wish to study two-mode nonclassical states and their potential to violate Bell-type inequalities then we can 
choose the other two modes to be in the vacuum state or in general in a classical state, then perform a suitable C/(4) 
transformation and proceed with the analysis. If we observe any violation, it can then be attributed to the initial 
nonclassical two-mode state. This is very similar to the detection of squeezing, where the mixing of the squeezed 
signal with high intensity coherent light is required in order to measure squeezing. 

Consider a 4-mode state with a general centered Gaussian distribution ||l^ 



= ( 91 92 <73 Qi Pi P2 P3 Pi ) 

G = real symmetric positive definite 8x8 matrix. 



(3.14a) 
(3.14b) 



Here g's and p's are quadrature components corresponding to a' and a^'s (gi = "^(oj + ai),?'i — "^('^1 ~ "^i) etc.). 
For W{^) to represent a quantum mechanical state and to be a Gaussian Wigner distribution, G has to satisfy in 
addition the condition 



f3- 



G-i > 

04x4 l4x4 
— I4x4 04x4 



(3.15) 



This is an expression of the uncertainty relations between the canonically conjugate q's and p's. The matrix V — 
is the variance or the noise matrix. For a given state, if the least eigenvalue £<:(V) of this matrix is less than ^ 
the state is squeezed and hence nonclassical MM : 



1 /^-i 
^ then 



S<{V) < ^ Squeezing 



(3.16) 



We now proceed to study the possibility of violation of the Bell type inequality ( |2.18| ) for Gaussian states and its 
possible correlation with squeezing. To calculate the coincidence count rates we need to express the J7(4) matrix 
corres pon ding to the rotation to the basis defined by the polarisers, as a matrix of Sp{8, R); we construct it directly 
from §J) 



C/(01,02) 



/ cos 01 
sin 61 






V 



— sin 61 
cos 01 












cos 62 
sin 62 














— sin 02 
cos 6*2 












cos( 
sin^ 










— sin 6*1 
cos 01 












cos 02 

sin 02 










— sin 02 
cos 02 



(3.17) 
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In order to calculate the coincidence count rates, we first compute the overlaps of the 4-mode Gaussian state with 
appropriate vacuum states. The Wigner function for a single-mode vacuum state is given by 



Wo{q,p) = -exp(-q2 

TT 

The relevant overlaps can then be written as]^ 

Tr(p|0)e, e^m^2n ( 0) W^o(9i,Pi)de 



= 2^Dei{UT{ei,Q)GU{di,Q) + en + 655)-! 



= 2y^Dct(i7^(0, 02)GC/(0, 62) + 633 + 677)-! 



= 4^Det(C/^(0i, 02)Gf/(0i, 02) + en + 633 + 655 + 677)^1 
Tr(p|00)kk(00|) = (2^)2 / M/(OW'o(<?i,pi)W^o(g2,P2)de 



= 4v/Dct(G + en + 622 + 655 + eee)-^ 
TV(p|00)k'k'(00|) = (2^)2 / W^(e)Wo(<?3,P3)W^o(g4,P4)de 



= 4VDet(G + 633 + 644 + 677 + ess)"^ 
Tr(p|00)kk(00||00)k'k'(00|) = {2Tif f (0 Wo(<zi,Pi)W^o(<Z2,P2)W^o(g3,P3)W^o(g4,P4K 



= 16v/Det(G + /8x8)-i 
Tr(p|0)9, e,(0||00)k' k'(00|) = {27rf J W {U{9i,0) OWo{qi,Pi)Woiq3,P3)Wo{q^,Pi)d^ 

= 8^Det(i7^(0i, O)Gf/(0i, 0) + /sxs - 622 - ege)"! 

Tr(p|0)e, e.(0||00)kk(00|) - (2^)3 ^ (t/(0, ^2) W^o(g3,P3)W^o(gi,Pi)Wo(<?2,P2K 

- 8y^Det(i7^(0, 02)GC/(0, 62) + /gxs - 644 - ess)-^ 



(3.18) 



(3.19) 



Using these overlap integrals, and the definitions ( 2.10| ) we can immediately compute the quantum mechanical coin- 
cidence count rates 

^(^i,^2)|r'' = 1 - Tr(p|0)e, gAO\) - Tr{p\0)g, ,,(0|) + Tr(p|0)e, 9i(0||0)e, g,{0\) 
m, )|r'' = 1 - Tr(p|0)e, e,(0|) - Tr(p|00)k'k'(00|) + Tr(p|0)e, (0| |00)k' k' (00|) 
Pi ,^2)fr'' = 1 - Tr(p|00)kk(00|) - Tr(p|0)e, e.(0|) + Tr(p|00)k k(00| |0)e, (,.(0|) 
^( , )|r''' = 1 - Tr(p|00)kk(00|) - Tr(p|00)k' k'(00|) + Tr(p|00)k k(00| |00)k' k'(00|) 

(3.20) 



The equations ( |3.19 ), can now be used to analyse any Gaussian Wigner state to see if it violates the inequality ( ^.18 ). 
Having set up the general formalism for the family of centered Gaussian Wigner states, we now consider examples of 
G which lead to violation of the inequality (2.18). Consider 



^Cij is an 8 X 8 matrix with {eij)ij — 1 and all other elements zero; /gxs is eight dimensional identity matrix 
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G = U-^S'^GoSU 

Gq = k/8x8, < K < 1. 



tanh— , B 

2' ' 



TlLO 

kf 



(3.21) 



Here k = I implies zero temperature and k < 1 corresponds to some finite temperature, S" is a 4-mode squeezing 
symplectic transformation, which is a Sp(8, 5R) matrix, and [/ is a passive symplectic J7(4) transformation whose role 
is to produce entanglement. As an example, we start with a state in which the modes ai and 04 are squeezed by 
equal and opposite amounts u and the modes 02 and 03 are squeezed by equal and opposite amounts v, and the 
entanglement is "maximum" . This corresponds to the choices]^ 



S 



/ 




























































































e" 






























e" 




























































e" 







v 

























J 



,u 



/ 


1 


1 


1 


1 















-1 


1 


-1 


1 
















-1 


-1 


1 


1 
















1 


-1 


-1 


1 




























1 


1 


1 


1 
















-1 


1 


-1 


1 
















-1 


-1 


1 


1 


v 














1 


-1 


-1 


1/ 



(3.22) 

For this particular class of centered Gaussian Wigner states the function (dependences on u and v are left implicit) 



yauss 

7qm 



(3.23) 



can be calculated. 



^This S matrix is constructed from two two-mode squeezing transformations studied in detail in [g^, the first being the one 
which squeezes modes 1 and 4 by equal and opposite amounts it, whereas the second one squeezes the modes 2 and 3 by equal 
and opposite amounts v 
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0.25 
0.2 

5 0-15 
CO 0.1 
0.05 

CO 

^ 



-0.05 




(c) 



-8 -6 -4 -2 2 4 6 8 
U=V 



CO 
CO 

oo 



0.1 



0.05 



■0.05 



-0.1 



(d) 



-8 -6 -4 -2 2 4 6 8 
U= -V 



FIG. 2. Violation of Bell type inequality for states with centered Gaussian Wigner distributions representing a 4-mode 
squeezed vacuum, (a) v = i.e., Modes 1 and 4 are squeezed by equal and opposite amounts w (b) w = i.e., Modes 2 and 3 
are squeezed by equal and opposite amounts {c) u = v i.e., all modes are squeezed, (d) u = —v i.e., all modes are squeezed. 
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4 6 8 



FIG. 3. Violation of Bell type inequality for states with centered Gaussian Wigner distributions representing 4-mode squeezed 
thermal states with the same temperature for all the modes. The choice of parameters for figures (a), (b), (c) and (d) is the 
same as that for Figure 2. 

Though one ought to search over all values of the angles 61,62, 6[ and 62 to look for possible violations of the 
inequality ( 2.1^ ), motivated by the choice of angles for the two-photon state, we restrict ourselves to that same choice, 
and plot the function /(-g-, ^, 3^, 0) for various combinations of the squeeze parameters in Figures 2 and 3. In Figure 2 
we have considered the case k = 1 i.e., the state under consideration is the squeezed vacuum; whereas for Figure 3 we 
have AC — 0.75 which corresponds to a squeezed thermal state. For the case of the squeezed vacuum, a clear violation 
is demonstrated in figures 2(b) and 2(c), for a considerable range of squeeze parameters. No violation is seen for the 
parameters chosen in Figure 2(d), and Figure 2(a) shows very little violation. However, for these parameter values 
for squeezing the violation might occur for some other values of angles 61,62, 6[ and 62. Figure 3 has similar features 
as Figure 2 but we see as expected that the amount of violation has diminished when the initial state is thermal 
instead of vacuum; as a consequence the small amount of violation which was present in Figure 2(a) has disappeared 
in Figure 3(a). As a further study, in the context of Gaussian Wigner states, it will be interesting to see the effect of 
phase space displacement on the violation of these inequalities; this will be taken up elsewhere. 
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IV. CONCLUDING REMARKS 



We have developed the machinery for analysing the violation of Bell type inequalities for a general state of the 
4- mode radiation field in a setup of the type described in Fig 1. A classical state in the quantum optical sense 
always obeys these inequalities while a nonclassical state may violate them, possibly after a ?7(4) transformation. 
Starting with a general nonclassical state, we subject it to a general unitary evolution corresponding to passive 
canonical transformations U (4) before we look for the violation of Bell type inequalities. It may turn out that a given 
nonclassical state does not violate Bell type inequalities but some U{A) variant of it does. 

A pure quantum mechanical state of a composite system is said to be entangled if we are not able to express it 
as a product of two factors, one belonging to each subsystem. Such states have nontrivial quantum correlations and 
can lead to the violation of suitable Bell's inequalities. The simplest quantum optical system for which the notion of 
entanglement can be introduced is the two-mode field. The group of passive canonical transformations in this case is 
C/(2); its elements, though incapable of producing or destroying nonclassicality are capable of entangling(disentangling) 
originally unentangled(entangled) states. As an simple example we choose the nonclassical state |1)|1) which is not 



it: , 



— {a\a2 + a\ai) -i 
entangled; by a simple U{2) transformation e 4 it becomes the entangled state ^^(|2)|0) + |0)|2)). 

v2 

However, coherent states are not entanglable in this way! If we start with a two-mode coherent state |zi,Z2), it 
is clearly not entangled i.e., the state can be written as a product with one factor belonging to one mode and the 
other to the other mode. Under a U{2) transformation this property is maintained. Classical states are statistical 
mixtures of coherent states and under U (2) transform again to such mixtures of coherent states. Such a mixture can 
definitely have correlations which are purely classical, but it cannot have truly quantum mechanical entanglement. 
Thus classical states are to be regarded as not entangled, and they remain so under passive U{2) transformations. 
However this is in general not true for a nonclassical nonentangled state which may get entangled under a suitable 
U (2) transformation. It is a straightforward matter to generalise the above statements to n mode systems where 
the group of passive canonical transformations is U{n). There are several ways to quantify entanglement; for pure 
states it is unambiguous: if the reduced density matrices for the subsystems involved are also pure states then the 
original state is not entangled. On the other hand, if in the process of partial tracing, some information is lost then 
the original state is entangled. The generalisation to mixed states is nontrivial but is conceptually simple; we have to 
separate classical correlations from the quantum mechanical ones and this may not always be easy to do. However, 
as we saw for the case of classical states, we can sometimes easily conclude that a given state is nonentangled. 

The above conclusions have an interesting bearing on the work on violation of Bell-type inequalities with beams 
originating from independent sources |Q ||l^. These experiments take two beams from two independent sources, 
pass them through some passive optical elements and show that the Bell-type inequalities are violated. The first 
conclusion we can draw from our analysis is that it must be the quantum optical nonclassicality of one of the beams 
in this experiment which has been converted into entanglement by the C/(4) transformation and hence led to the 
violation. Secondly, if the original beams were quantum optically classical, no matter what one does, no violation will 
be seen. 

In our analysis, we have not distinguished between strengths of coincidences. The coincidence counter registers a 
count when simultaneously each detector detects one or more photons. This is the reason why we chose the operators 
A's to have eigen values and 1. In this sense, the measurements involved here are not refined. It would be interesting 
to further generalise the analysis by considering somewhat refined measurements where to some extent coincidences 
are distinguished on the basis of their strengths. However, the relevant operators in this context may be unbounded; 
and it is well known that the formulation of Bell type inequalities for such operators, though desirable, is nontrivial. 

We have compared quantum optical nonclassicality with violation of Bell's inequalities. When a state is nonclassical 
in the quantum optical sense, it does not allow a classical description based on an ensemble of solutions of Maxwell's 
equations, which is a very specific classical theory. On the other hand violation by a state of a Bell type inequality 
rules out any possibility of describing it by any general local "classical" hidden variable theories. Therefore, it is 
understandable that quantum optical nonclassicality is a necessary but not a sufficient condition for the violation of 
Bell's inequalities. This disparity is partially compensated for by the freedom to perform passive canonical transfor- 
mations on a nonclassical state before looking for violation of Bell's inequalities though it is not obvious whether this 
freedom completely removes this discrepancy. On the other hand, if a state obeys Bell's inequalities, it may still not 
allow a "classical" description. Therefore, we need a complete set of Bell's inequalities capturing the full content of 
the locality assumption. These and related aspects will be explored elsewhere. 
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